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Abstract 

We study the effects of scalar-pseudoscalar mixing induced from quantum loop on e'/ e in 
SUSY models with tan/? ~ nit/nit,. We find that even the non-universal soft A d term and 
Yukawa matrix, Y d , are hermitian, the predicted value of |e'/e| can be consistent with the 
measured results of NA48 and KTeV. And also the EDMs are compatible with experimental 
bounds. 



Since its discovery in the neutral kaon decays in 1964 [I], the origin of CP violation 
(CPV) still puzzles the physicist. Although recently another time-dependent CP asymmetry 
(CPA), sin 20 1; in the decay of B -> J/^K S is observed by BARBAR [g and BELLE f|, our 
understanding of CPV is not much better than before. What we are certain at present is only 
that it is necessary to exist unrotatable phase and it is believable that such kind of phase is 
associated with weak interactions, usually called weak phase. In the case of standard model 
(SM), the unique source of CPV is from the Cabibbo-Kobayashi-Maskawa (CKM) matrix 
Q induced from the three-generation quark mixings and described by the three angles a, (3 
and 7 or 02, 0i and 03. 

Even though the SM prediction on the indirect CP violating parameter e in the kaon 
system can be fitted well with current experimental data, due to the large uncertainties from 
hadronic matrix elements, so far it has not been settled yet whether the result in the SM 
can explain the observed value of the direct CP violating parameter e' measured by NA48 
H and KTeV |5|]. And also, because of the unitarity in the CKM matrix, the predicted 
electric dipole moments (EDMs) of neutron and lepton are quite small and unreachable 
experimentally. Moreover, the requirement of Higgs boson with the mass being less than 60 
GeV to solve the problem of baryogenesis is ruled out by the LEP experiment. Therefore, 
it becomes important to search the possibility of existing other CP violating sources for 
explaining all CP phenomena. 

One of reliable models in the extension of SM is supersymmetric (SUSY) model. SUSY 
theories not only supply an elegant mechanism for the breaking of the electroweak symmetry 
and a solution to the hierarchy problem, but also guarantee the unification of gauge couplings 
at GUTs scale J7|. In addition, SUSY possesses abundant flavor structures, such as upper 
and down type squark mixing matrices, and CP violating phases, which are arisen from the 
trilinear and bilinear SUSY soft breaking A and B terms, the \i parameter for the scalar 
mixing etc.. Unfortunately, one can check easily that those phases are severely bounded by 
electric dipole moments (EDMs) || so that the effects on e and e' cannot enough explain the 
current experimental values. In order to handle the small CP phase problem, it has been 
suggested to use the non- universal soft A terms instead of universal ones ||. Furthermore, 
to evade any fine-tuning on the specific phases which contribute to EDM of neutron, it 



is proposed that SUSY soft-breaking A g and Yukawa, Y q , matrices are hermitian |I0], |TT 



Consequently, the CP phases of 0(1) can exist naturally even considering the contributions of 
EDM. And also, it implies that the CPA in hyperon decays could reach the value of O(10 -4 ) 
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p~2|| proposed by the experiment E871 at Fermilab |[L3|| . However, the effect on e', dominated 



by gluon-penguin with gluino as the internal particle in the one-loop, will be suppressed due 
to (S( 2 )lr — (^12) rl ( the definition is shown below). It is interesting to ask whether e', 
based on the gluino mechanism, can be satisfied in the framework of hermitian A d and Y d 
matrices by considering other effects which are available after including the constraints from 
experimental measurements. 



According to the analysis in |]14 |, we find that because of the enhancement of large A q 1 
fi and tan f3, the coupling iV° — qi — qn, with N° and q being scalar (or pseudoscalar) and 
the squark of corresponding to the g-quark respectively, will be enhanced. Thus, it is easy 
to conjecture that the vanished e' problem in hermitian case could be solved if we consider 
the induced interactions Cf^oN f Ri f L j in which Cf^o is the effective coupling and / could 
be upper or down type quarks and the indicies % and j denote the possible flavour. From the 
induced effective vertices, if Cjn° is complex, we see clearly that e will also contribute. If 
so, we will face the strict constraints from it. In our following analysis, the induced effective 
coupling will be taken as real. But, the CPV is generated by the scalar-pseudoscalar mixing. 

It is known that if the scalar-pseudoscalar mixing comes from spontaneous CPV, such 



as realized by radiative corrections in the MSSM |jl5|j , the predicted pseudoscalar mass, m^, 
is far below the current experimental limit and excluded |]IE|| . However, it is found that if 
CP is broken at the tree level by the SUSY soft breaking sector, the large mixing between 
CP-even and CP-odd boson can be obtained through radiative corrections without the limit 



of small mass on the pseudoscalar boson [17|. As also been shown in [18|, the radiative CP 



effects will modify the couplings of Z-boson to Higgs boson such that the mass of lightest 
Higgs boson with 60 — 70 GeV can escape from the bound of LEP2. With above conclusions, 
in this paper, we will show their implications on e'. 

We start by writing the effective interactions between the relevant neutral Higgs bosons 
or gluino and squarks as 

£eff = £<N~f L f R + 



2M W 

-iA°&((Ml) LR Z$- ^* mq )q R 

-V2g s [qP R g a T a q L - qP L g a T a q R ] + h.c. (1) 

where Pl(r) = (l^Ts)^, h° stands for the lightest scalar particle while A is for pseudoscalar 
boson. The bold q and q denote three generaton quarks and the corresponding squarks, 



respectively. The generators of SU(3) C are normalized by tr(T a T b ) = l/2S ab . = 
cos a/ sin/5 (—sin a/ cos 0) and 0)3 = — sin a/ sin/5(coso;/ cos/5) in which angle a describes 
the mixing between two CP-even Higgs particles. For simplicity, we only concentrate on the 
contributions from the lightest CP-even Higgs. Due to the mass suppression, we expect that 
the effects of heavier CP-even boson are smaller than those of light one. One the other hand, 
the squared squark mass matrices responsible for flavor change are described by 

2 _ / (m?) LL (m?) M \ 



m 2 ) LL = (Ml) LL + m 2 -M 2 zC os2pC q L i, 



q 

(m 2 g) LR = (M 2 ) LR - fi*Z q pm q , 

(m 2 ) RR = (M*) RR + m 2 q + M 2 z cos2f3C R i, (3) 

where we have adopted the so-called super-CKM basis that the quarks have been the mass 
eigenstates so that m q is the diagonalized quark mass matrix, q and q stand for quark and 
its superpartner. They could be upper or down type quark. Z^' = cot/5(tan/3) and 

C q L = T^Q q sm 2 6 w , 

C q R = Q q sm 2 9 w (4) 

with T 3 and Q q being the z-component of isospin SU{2)l for the squark q and its charge, 
respectively. 1 denotes the 3x3 unit matrix. The definition angle (3 is followed by tan (3 = 
v u /vd with v u and Vd being the vacuum expectation values (VEVs) of Higgs fields and 
$ d responsible for the masses of upper type quarks and down type quarks, respectively, 
is the mixing effects of $ u and $ d . (M-^)ll(rr) stand for the soft breaking masses for the 
corresponding squarks and (M?)^ describe the trilinear soft breaking couplings and are 
written as 

(M?) M = ^V qL A^ R , (5) 

where V qL{R) transform the left(right)-handed quarks from weak eigenstates to mass eigen- 
states and Ajj = Y q jA% with being the Yukawa matrix. 



For convenience, we adopt the so-called mass- insertion approximation method [IS 
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which the basis for squark is chosen such that the gluino- squark- quark vertices involving 
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quarks are flavour diagonal (super-CKM basis) instead of diagonalizing the squark mass 
matrix itself. Hence, the squared squrak mass matrices are regarded as effective couplings. If 
necessary, we can insert the proper effective couplings in the propagator of squark. According 



n 9 

\q Ai q Bj / 



N° 



(a) 



g 

q Ai \, ~ ~ t /q B j 



N c 



(b) 



iBj 



Figure 1: Feynman diagrams for one- loop induced couplings of neutral Higgs to quarks: 
N° could be scalar or pseudoscalar boson, q could be upper or down type quarks, q is the 
corresponding superpartner, and the indices i,j,k stand for the possible flavour of quarks. 
AO and £0 denote the chirality and = R(L) while = L(R). 



to Eqs. ([I]) and (0), the interactions for N°d B s A and N°q B q A , illustrated in Figure [I], can 
be derived as 



C = (AV2G F ) l ^[(N^) BA d B s A N° + (N^) BA d B d A N 
+ (N^) BA u B u A N° 



(6) 



where A and B denote the chiralities and they are always opposite to each other, (N g ) BA - 
(H Q ) BA ((A q ) BA ) if N° is scalar (pseudoscalar) Higgs and their expressions are written as 



(N, 



ij 



IBA 



Arc 



n 



Bij 



+ ~( S kj)B'Bn B , kk (5* k ) AA/ I 2 (x Q 



(7) 



0$ 



)AA' 



AA> 



with x q = m 2 g/m 2 q and m g and irtq being the average masses of gluino and q type squark. 
The first term in Eq. (|7]) is from the lowest order contributions, illustrated in Figure |](a), 
while the second term is generated by double mass-insertion, shown in Figure [l|(b). The 
chirality A' in Eq. @ can be L or R, but chirality B' is opposite to A' so that beside 
i^fj) lr(rl) j ll(rr) will also contribute. fi\ = h q A {d q A ) are related to the couplings of 
scalar (pseudoscalar) boson to the squark q in which the expressions are 



hi 



- q )LR q n m q 

^1 H ^2) 



rrig 



and n q R = fl 9 ^ . 



al = J^^Z q +i^ 
nig nig 



f . 1 x\u{x 
Ux) = + 



1 — X (1 — X 



2 ■ 



T , , 2 + 3x — (6 — x)x 2 3x ln(x 

= — — — \- 



2(1 -x) A (1-x) 4 

We note that the relevant effects in the first term of Eq. (|7|) are h d L12 , h d Lll , and h u LXV Via 
their definitions, we have 

d ™ / rd \ ,~.d , A 1 * 77 ^! ~ d 



^11 = rn i (5f 1 ) LR Pf + 



to. 



m d 



2- 



with m u i and m^i being the masses of u-quark and d-quark, respectively. According to the 
analysis in |[20|| , the constraint on |/m((5 d 2 )LR| from e' is of order of 10~ 5 . With the assumption 
of the phase of 0(1), we expect that the |(<^ 2 )l-r| has the similar order of magnitude so that 
its contribution is negligible. Although it is not necessary, to simplify our analysis, we adopt 
that the Yukawa matrices have the structure > Y^ d ^ in which Y^ d ' is any entry except 
Y 3 3^ d ' ) and each of them is proportional to ^Jm^m^j m 3 with m 3 being the top or bottom quark 
mass for the corresponding Yukawa matrix f2~lf . As a result, {&\\)lr will be related to the 
first two generation quark masses or be suppressed by flavour mixing elements as defined 
in Eq. (|5|). In this paper, we also neglect their contributions. The similar situation is also 
applied to a q L . According to the above assumption, we see that the dominant effect in second 
term of Eq. (|7]) is from the component n| 3 which are expressed by 

~ h L33 = mq(5^ 3 ) LR ul + — — ^>g, (8) 

TTlq 

«L33 = -imq(8l 3 ) LR Z q p + i—^-. (9) 

Wig 

It is worth mentioning that due to hf^^ and a^ 33 associated with 1 / cos /3 and tan j3 respec- 
tively, if considering large tan/? case, there is a large enhancement. 

In terms of Eq. @, we can immediately obtain the effective operators for IAS*! = 2 and 
|A5| = 1 as 

(H 12 )ba(H 12 )b'A' 



C\as\=2 = -AV2G F Yl { 



2 

A,A'=L,R " m h 



AS|=1 



xA 5P J d B SA d B is A i 
-4y/2G F Yl \ 



(10) 



(H 12 )ba(Aii)b'A' 



A,A'=L,R 

(M2)ba(Hi 1 )b'a> 



mi 



nit 



AspQbQA d>B'SA' 



where q = u,d and A$p = Mg P /m A describes the scalar-pseudoscalar mixing effect and its 
qualitative dependence can be described by [[IT], [H| 

m\ Im{A?ii) ( \A l \ 2 



\ 2 

l SP 



2Re(A t fi) 

^ 327T* MI USY [ L > MI USY > tan Wh^' ^susy , 



(12) 



In terms of the results in one can find that one loop radiative effects M$p could be 
of order of few hundred GeV, that is the CP mixing factor A$p could be 0(1). In order 
to study the contributions to \AS\ = 2 and 1 decays, such as Am K , e and e', the relevant 
hadronic matrix elements are estimated by vacuum saturation method and written as 



(K°\d R s L d R s L \K°) 

5 / m K f 

24^m s + mJ ™ 
(K°\d L s R d R s L \K°) 



r 1 




/ mjc \2- 


-24 " 




"> m s + md ■ 



ttik/k 



(13) 



for AS = 2 decays |22| and 

V 1 = (n-\d l5 u\0)(n + \us\K ) = -f v B 2 (l + 2 l ^f 

V 2 ~- 



'A 2 

x 

mh 



7T 7T 



|gg|0>(0|J 7 5^ > = -f*B*(l + 2jf) 



(14) 



with B = m\j '(m s + md) and A x « 1. GeV for | AS 1 ) = 1 |23|]. What we are concerned is only 
the main effects on the FCNC decays, thus, due to mh < m A , in the following calculations, 
we only concentrate on the contributions of in Eqs. (|I~0|). 
It is known that the observable for e'/e is expressed as 

Re (i) a 



■(Q 1 ImA 2 - ImAo) 



(15) 



/ V2ReA \e 

where A and A 2 denote the two-pion final states of — > tttt in isospin 1 = and 1 = 2, 
respectively, Q = ReA 2 /ReA w 1/22 and ReA w 2.7 x 1(T 7 GeV. From Eqs. (JTT|) and flU), 
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we know that the new effects on ImA 2 and ImA have similar value in magnitude so that due 
to the suppressed factor Q in isospin 1 = 0, the final state with / = 2 is dominant. Because 
the CP violating phases are arisen from soft SUSY breaking terms, in our considering case 
the ll(rr) are real. On the other hand, although (^ 3 )l_r(_rl) (i = 1, 2) could be complex, 
however, charge and color breaking (CCB) minima and the potential unbounded from below 
(UFB) will give strict constraints so that their contributions are small and negligible. In 
order to avoid the constraints from e, we impose that Arg(fi) = 0, A d and Y d are hermitian 
matrices so that {5 33 )lr{rl) is real, and (S^ll ~ i^ij) RR ■ As a result, h\ 33 and d^ 33 are real 
and purely imaginary, respectively. And then we get the identities {Hf 2 ) LR = (Hf 2 ) RL and 
(Af 2 ) LR = -(Af 2 ) RL such that the e relating effect, such as (Hf 2 ) RL (Af 2 ) RL + (Hf 2 ) LR (Af 2 ) LR 
is vanished. Altogether, Eq. (|10|) can only contribute to Am^. Moreover, from Eq. fljjp, 
we clearly see that the M| P is related to 7m(A*/x), therefore, (5 33 )lr in Eq. (^) should be 
complex. Nevertheless, in terms of Eq. (|TTD and due to purely imaginary (Af 2 ) R A, it is 
obvious that only real part has the contribution. 



Combining the results of Eqs. (ITO) — (O), we get 



Am K = 2Re(K \n lAS \=2\K ) 



16v^GF ^ffi (£i + K 2 ), (16) 



m 



,e'\ 2V2G F , 1 v 

Re \e) I=2 ~ 3mlReA \e\^ ~ 2Nj Vl 



x 



^w^nW + (i^fltWOflLlAsp (17) 



where iV c = 3 is the color number, for simplicity we have used (Hfj) R L ~ (H^-)xi? and 
(Alj) R L ~ — (A^-)lr for q = u, d. Because the constraint on (5^ 3 )ll is stricter than 
that on (5± 3 )ll, in Eq- (H3)> we °nly show the contributions of (H^ab an d (^u)ab- 
To escape the constraint from Atjik directly, we set the A b , \i and angle a satisfy with 
cot a « m'~(5 33 ) R l / firrib so that (Hf 2 ) R i ~ 0. In order to obtain the measured value of 
e', the values of relevant parameters are taken as « = m|, ~ x u = = 0.3, 
(<%) M ~ 0.1m g ~/500GeV, (<%) z£ ~ 0.3m f /500GeV 0, («^ 3 ) LL ~ 0.4 (m 9 ~/500G e y) 2 , 
~ 4.5 x 10- 2 K~/500GW) |2|, = \A b \/m q ~ 2, /i/m g ~ ~ 1, (5^ d) ) RL « 

A^'rub/m^, rriq « 800 GeV and tan/? ~ m t /mb. As mentioned before, the magnitude of 
squared mass arisen from radiative corrections for the mixing between CP-even and CP-odd 
boson could be order of (lOOGeV) 2 , i.e., A S p could be order of unity. As a consequence, from 
Eq. fll7l) and above taken values, the predicted direct CP violating parameter for — > nn 



is given by \Re(e'/e)\ » 1.83 x 1(T 3 with m h = 120 GeV and A SP « 0.30. The result is 
consistent with (15.3 ±3.6) x 10~ 4 and (20.7±2.8) x 10~ 4 given by NA48 § and KTeV §, 
respectively. 

Finally, we give the estimation on the electric dipole moments (EDMs) of electron and 
neutron. In our present considering case, the scalar-pseudoscalar mixing only depends on 
the complex (S 33 )lr- That is, even the mixing between CP-even and CP-odd boson is small, 
it still can introduce CP violating effects, such as EDMs of neutron and lepton. According 



to the results of [26], the neutral Higgs can contribute to EDMs of neutron and lepton via 
two-loop topologies. Due to Asp being less than unity, the main effects should be from 



pseudoscalar exchange. Hence, following the results of |26| , the EDM of fermion can be 
written as 

fey = Ql N £ 7 _ (is) 

\ p. I J 327r 3 m\ L \ m 2 A J ^m\^ 



where t\ and ti are the mass eignestates of stop-quark, Q f denotes the corresponding fermion 



charge and £ 4 = Z^[im t Im(5 33 )LR/ sin/5 cosfiv 2 with v = yjv% + vj, the definition of func- 
tion F can be found in |2"6| . By using naive valence quark mode and taking g s (A) = 47r/v^6, 
a s {M z ) = 0.12, m u = 7 MeV, m d = 10 MeV, m A = 400 GeV, the predicted EDMs of 
electron and neutron with QCD renormaliztion effects, described by 

are around 1.01 x 10~ 27 cm and 3.62 x 10 -27 cm, respectively. Both are satisfied the current 
experimental limits given as d e /e < 4.3 x 10~ 27 cm [[Z7| and d^fe < 6.3 x 10~ 26 cm |28|| . 

In summary, we have studied the effects of scalar-pseudoscalar mixing on e'/e in SUSY 
models with tan/? ~ mt/m^. We find that if the non-universal soft A d term and Yukawa 
matrix, Y d , are hermitian, and (<5^-)ll ~ (Sfj)RR, the predicted value of e'/e is consistent 
with the measured results of NA48 and KTeV. And also the EDMs are compatible with 
experimental bounds. 
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